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Abstract 

We study three regimes of the Mott transitions characterized by classical, marginally quantum 
and quantum. In the classical regime, the quantum degeneracy temperature is lower than the 
critical temperature of the Mott transition, T c , below which the first-order transition occurs. The 
quantum regime describes the T c = boundary of the continuous transition. The marginal quan- 
tum region appears sandwiched by these two regimes. The classical transition is described by 
the Ising universality class. However, the Ginzburg-Landau- Wilson scheme breaks down when the 
quantum effects dominate. The marginal quantum critical region is categorized to a new univer- 
sality class, where the order parameter exponent (3, the susceptibility exponent 7 and the field 
exponent 5 are given by j3 = d/2, 7 = 2 — d/2 and 6 = 4/d, respectively, with d being the spatial 
dimensionality. It is shown that the transition is always at the upper critical dimension irrespective 
of the spatial dimensions. Therefore, the mean-field exponents and the hyperscaling description 
become both valid at any dimension. The obtained universality classes agree with the recent ex- 
perimental results on the Mott criticality in organic conductors such as k-(ET)2Cu[N(CN)2]C1 and 
transition metal compounds such as V2O3. The marginal quantum criticality is characterized by 
the critically enhanced electron-density fluctuations at small wavenumber. The characteristic en- 
ergy scale of the density fluctuation extends to the order of the Mott gap in contrast to the spin 
and orbital fluctuation scales and causes various unusual properties. The mode coupling theory 
shows that the marginal quantum criticality further generates non-Fermi-liquid properties in the 
metallic side. The effects of the long-range Coulomb force in the filling-control Mott transition are 
also discussed. A mechanism of high temperature superconductivity emerges from the density fluc- 
tuations at small wavenumber inherent in the marginal quantum criticality of the Mott transition. 
The mode coupling theory combined with the Eliashberg equation predicts the superconductivity 
of the d x 2_ y 2 symmetry with the transition temperature of the correct order of magnitude for the 
realistic parameters for the cuprate superconductors. Experimental results on the electron dif- 
ferentiations established in the angle-resolved photoemission experiments are favorably compared 
with the present prediction. The tendency for the spatial inhomogeneity is a natural consequence 
of this criticality. 

PACS numbers: 
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I. INTRODUCTION 



Mott transition belongs to one of the metal-insulator transitions ubiquitous in various 
compounds 1 . Physical properties of Mott transition and its nature are a long-standing 
subject of research with many controversial issues. The problem has first been postulated in 
1937 by Peierls and Mott 2 . Although the band theory of metals and insulators established 
soon after the foundation of quantum mechanics is quite successful, it was pointed out by 
de Boer and Verway^ that insulating behaviors of many transition metal compounds as 
NiO cannot be explained by a simple band picture, because the bands are partially filled. 
Peierls pointed out a crucial role of electron correlations as the mechanism of the insulating 
behavior. Mott developed this idea and introduced a concept which we nowadays call the 
Mott insulatoi^. Since then, it has been recognized for a long time that the Mott insulator 
stabilized by the electron-electron Coulomb interaction and metallic states near it provide 
us with fruitful physics with various novel concepts. This has become very popular after 
the discovery of the high-T c cuprate superconductors, which was indeed discovered in doped 
Mott insulators^. However, relationship of this fruitful outcome to the nature of the Mott 
transition itself has not been fully clarified. 

The Mott transition can be realized basically by two routes: Bandwidth-control and 
filling-control. In the first route, the bandwidth is controlled relative to the amplitude of the 
local electron-electron interaction by keeping the electron density fixed at a commensurate 
value (namely, the electron density per unit cell is kept at an integer value). This route may 
be experimentally realized by applying pressure or by substituting elements with a different 
ionic radius and the same valence. In this route, the overlap of the wavefunctions between 
neighboring electronic atomic orbitals forming conduction electron bands is controlled. In 
the second route, electron filling is changed from the Mott insulator. This is typically realized 
by substituting with elements, which have a different valence from the substituted elements 
in the reservoir structure of the Mott insulator. 

Originally, Mott has considered only the first route, the bandwidth-control transition. 
In this category, there exist several typical Mott transitions including those observed in 
V2O3, i?Ni03 and k-ET type organic compounds, where R represents a rare earth element. 
Many of them show first-order transitions between the Mott insulator and metals at low 
temperatures^. In many cases, some magnetic and/or orbital order exist at low temperature 
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FIG. 1: Schematic and typical phase diagram of metal and Mott insulator in the parameter space 
of temperature T and the interaction strength scaled by the transfer, U/t. The solid curve shows 
the first-order Mott transition with the critical end point illustrated by the solid circle. In many 
cases, the antiferromagnetic (AF) phase appears at low temperatures as is illustrated by the broken 
phase boundary line. 

of the Mott insulator phase. Although sometimes concurrent magnetic transitions occur 
with the Mott transition at low temperatures, in the typical phase diagram illustrated in 
Fig. m the first-order Mott transitions occur even at high temperatures, where the mag- 
netic and orbital orders are not involved. This typical phase diagram clearly indicates that 
the Mott transition is inherently independent of the magnetic and orbital-ordering transi- 
tions. The statement that the mechanism of the Mott transition is primarily independent 
of the symmetry breaking of spins or orbitals is also corroborated by the existence of a 
quantum spin liquid phase in the Mott insulator recently found numerically^ as well as 
experimentally 9,10,11 . The Mott transition between a quantum spin-liquid and a metal does 
not accompany magnetic transitions. In fact, as we see in this paper, the Mott transition is 
not driven by the spin or orbital degrees of freedom but by the density degrees of freedom. 

Mott argued that the long-range part of the Coulomb interaction is necessary to reproduce 
the first-order transition 12 . The argument by Mott was the following: We can identify the 
Mott insulating phase as that where the two electrons sitting on the same atomic orbitals 
(we call it a doublon) and an empty site (we call it a holon) make a bound state. Metals 
are characterized by the phase where the bound state disappears. In terms of the binding 
of the doublon and holon, the binding energy is controlled by the screening of the attractive 
Coulomb interaction between the doublon and holon. Since the screening relies on other 
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doublons and holons moving in between, the screening becomes rapidly poor with decreasing 
density of free doublons and holons. Then at some threshold concentration of free doublons 
and holons, the screening becomes too weak to keep a doublon and a holon free and as a 
consequence of this feedback, they suddenly form a bound state, which leads to a first-order 
transition to a Mott insulating state. Another argument for the origin of the first-order 
transition emphasizes couplings to lattice distortion. When the metal is stabilized, the 
lattice constant diminishes further to gain the kinetic energy of electrons, which in general 
strongly favors the first-order transition through the electron-lattice coupling. 

However, recent detailed numerical studies on the Hubbard model on the square lattice 
overturned these speculations^^: Although the long-range Coulomb force and the coupling 
to the lattice may have some roles in the realistic Mott transition, even the Hubbard model 
with only the local onsite interaction without any coupling to lattice distortions is enough 
to reproduce the first-order Mott transition. The iV-site Hubbard model is defined as 

H = H t + Y^ H ui-vMN (1) 

i 

Kt = -J^tiMaCja + h.C.) (2) 

(ij) 

and 

n m = U{n^- l -){n ii -\). (3) 

Here, M = ^2 ia n ia /N and n ia = c\ a Ci a with the creation (annihilation) operator c\ a (ci a ) 
of an electron at the site i with the spin o. The chemical potential is /i and U is the 
onsite Coulomb repulsion. The phase diagram of the Hubbard model at zero temperature 
T = on the square lattice with the nearest-neighbor and next-nearest-neighbor transfers 
t and t', respectively, studied by the path-integral renormalization group metho d 14 ' 15 is 
shown in the plane of U/t and the chemical potential fi at zero temperature in Fig. |2J" 3 . 
The filling-control transition occurs across the edge of the boundary in Fig. El while the 
bandwidth-control transition is realized through the corner of the phase boundary at the 
bottom. The first-order transition through the bandwidth-control route is indicated by the 
jump of the averaged doublon density D = (n^n^) = Ylii 71 ^ 71 ^/^ as i n Fig. El for the 
case of the anisotropic triangular lattice^. When the first-order transition takes place at 
zero temperature, its jump decreases with raising temperatures and closes at the critical 
endpoint. In fact, the Hubbard model on this anisotropic triangular lattice is a relevant 
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FIG. 2: Phase diagram of Mott insulator (shaded area) and metal for two-dimensional Hubbard 
model on a square lattice with the nearest-neighbor transfer t and the second-neighbor transfer t' 
in the plane of the local interaction U and the chemical potential fi. The energy unit is t = 1 and 
t' is taken at 0.2. The bandwidth-control route is realized across the bottom corner of the Mott 
phase boundary along the gray dashed line, whereas the filling-control transition occurs across the 
solid lines. 

effective model for the k-ET type organic compound and the first-order Mott transition 
with the finite-temperature critical end point was observed experimentall y 16ll7lls . 

Another remarkable feature in Fig. |3] is that the jump in D decreases with increasing 
the next nearest neighbor transfer t'. This means that the so-called frustration effects re- 
duce the first-order jump and drives to more continuous type transitions. In fact, recent 
experimental studies on pyrochlore compounds* 9 appear to show a continuous Mott transi- 
tion by the bandwidth control with a signature of the Anderson localization in the vicinity 
of the boundary. We note that the pyrochlore lattice has the fully frustrated structure, 
where the antiferromagnetic order is severely suppressed. In addition, recent studies on k- 
(ET) 2 Cu2(CN) 3 suggest that the critical temperature of the Mott transition becomes lower 
for the more frustrated structure, namely for compounds with larger effective t' in the cor- 
responding theoretical model. — All of these consistently suggest that the order of the Mott 
transition may be systematically controlled from the first order to continuous. 

When the filling is controlled, the first-order transition appears as the phase separation. 
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FIG. 3: Double occupancy D at T = as a function of U/t for three choices of parameters, 
t'/t = 0.5, 0.8, and 1.0 for the anisotropic triangular lattice structure depicted in the inset with the 
nearest neighbor transfer t (solid bond) and the next-nearest-neighbor transfer t' (broken bond). 
Jumps of I? at U/t = 4.1,5.0 and 5.2 for t'/t = 0.5,0.8 and 1.0, respectively show the first-order 
Mott transitions. Curves are guide for the eyes 

Experimentally, the existence of the phase separation or the electronic inhomogeneity is 
a controversial issue as we discuss later. From theoretical side, numerical studies on the 
Hubbard model show the marginal result, where the phase separation does not occur while 
the charge susceptibility shows a critical divergence at the Mott transition . 21 ' 22 ' 23 In the 
terminology in this paper, we use the charge susceptibility and the density susceptibility as 
the same quantity d(M) /dfj, for the electrons with charge. 

Since the theoretical and experimental results suggest the controllability of the order 
of the Mott transition and its critical temperature, it is desired to understand the whole 
feature of the Mott transition from the universality classes of the finite-temperature critical 
point to the zero-temperature critical phenomena on the same grounds. We will show in 
this paper that three regimes of the Mott transition exist. One is the classical transition at 
a high temperature, which is described by the Ising universality class of the critical point 
accompanied by the first-order transition below the critical temperature. The second is 
the quantum transition, where the transition appears only at zero temperature, and the 
density susceptibility remains finite for the spatial dimension d > 2. The third regime 
is the marginally quantum one, which emerges at the crossing point of the classical and 
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quantum transitions. The marginally quantum regime is characterized by the diverging 
density susceptibility at small wavenumber for d > 2 at low temperatures. In the second 
and third regimes, the conventional scheme of the Ginzburg-Landau- Wilson theory does not 
apply. 

The new universality class at the marginally quantum transition has a deep consequence 
on the induced non- Fermi-liquid behavior in the metallic side, electron differentiation in the 
momentum space and the sensitivity toward electronic inhomogeneity. One of the most re- 
markable consequences is the superconductivity emerging from this marginal quantum Mott 
criticality. We show that the high temperature superconductivity of the d x -z_ y 2 symmetry is 
obtained under the realistic choice of the parameter values for the cuprate superconductors, 
where the density (charge) fluctuations at small wavenumber play the crucial role for the 
Cooper pairing. The energy scale of the fluctuation is characterized by the Mott gap, which 
can be by far larger than the energy scale of magnetic and orbital fluctuations for the filling- 
control transition. This solves many puzzling experimental results in strongly correlated 
electron systems particularly in transition metal compounds. 

A part of the discussions in this article is already give n 24 ' 25 . We summarize the pre- 
vious results and further extend the discussion on the quantum Mott criticality and its 
consequences in greater detail. In particular, detailed analyses on the breakdown of the 
Ginzburg-Landau- Wilson scheme for the quantum Mott transition are presented together 
with the scaling analysis. The validity of the mean-field exponents and the compatibility 
with the hyperscaling description are discussed in detail. Two-dimensional systems are espe- 
cially analyzed and are compared with the experimental results for the organic compounds 
and the cuprate superconductors. Basic finite temperature effects are also obtained and 
discussed in connection with the experimental results. An important issue for the filling- 
control transition is the effects of the long-range Coulomb interaction. We discuss how the 
present results are modified in the presence of the long-range interaction and also discuss 
the experimental relevance. The non-Fermi-liquid properties and the anisotropic Cooper 
pairing originating from the nonperturbative enhancement of the density susceptibility at a 
small wavenumber is an important subject we study in this article in detail. 

In Sec. II, we summarize the conventional Ginzburg-Landau- Wilson scheme of the Mott 
transition. In Sec. Ill, the nature of the quantum Mott transition is examined with emphasis 
on the breakdown of the Ginzburg-Landau- Wilson scheme. Section IV is devoted to the 
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resultant non-Fermi-liquid behavior expected in the metallic side of the critical region of the 
Mott transition. In Sec. V we discuss the mechanism of high-temperature superconductivity 
arising from the marginal quantum Mott criticality. Section VI concludes and summarizes 
the paper. 

II. CONVENTIONAL GINZBURG-LANDAU- WILSON SCHEME 

Recently, the critical endpoint of the first-order transition line of the Mott transition in 
Fig. [T] has been a subject of intensive studies. In case of V 2 3 , from the detailed study 
of the conductance, it has been suggested that the criticality of the transition follows the 
Ising-type universality class 26 . In an organic conductor of k-ET salt, the diverging electronic 
compressibility at the critical end point has been probed by the ultrasound velocity 16 . 

In prior to these experimental studies, several theoretical studies have focused on the 
nature of the transition. The Mott transition by itself does not change any symmetry. 
Therefore, from theoretical point of view, it has an analogy with the text-book gas-liquid 
transition, which is known to be equivalent to the ferromagnetic transition in the Ising 
model under magnetic fields. The first-order metal-insulator transition corresponds to the 
Ising transition between spin-up and down phases taking place with switching the direction 
of magnetic fields below the critical temperature. In fact, CastellaniS! has discussed the Ising 
nature of the first-order Mott transition by extending the Blume-Emery-Griffiths model^ for 
the phase separation of 3 He- 4 He mixture. It has further been considered in the dynamical 
mean-field theory by Kotliar et al . 29 ' 30 , where the low-energy part of the single-particle 
Green's function appears to follow Ginzburg-Landau scheme in accordance with the mean- 
field theory of the Ising model. 

In the Ginzburg-Landau- Wilson scheme for the Ising-type transition, the free energy may 
be expanded by the spatially dependent scalar order parameter X(r) integrated over space 
coordinate r as 

F = J dr[^a (T - T c )X{rf + \bX{r f - fxX(r)} (4) 

near the critical temperature T c with and h being positive constants. In the Ising model, 
X is indeed the order parameter, namely, the magnetization m. In the mapping to the gas- 
liquid transition, X is interpreted as the density of particles, n, measured from the critical 
density. At the critical temperature T c of the gas-liquid transition, the uniform density 
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susceptibility Xn = [d 2 F/dn 2 ]~ l diverges. 

When it is further extended to the mapping to the Mott transition, for the filling-control 
transition, X is identified indeed as the electron doping concentration X measured from the 
critical density at the critical point of the Mott transition^- This is a natural consequence, 
because, in the filling-control transition, the control parameter is the chemical potential, 
which is conjugate to the carrier density. 

When the bandwidth is controlled, the control parameter is U/t in the Hubbard model 
and the conjugate quantity to U is the doulon density D. Therefore, the order parameter 
X in this case is the doublon density D, which indeed jumps at the first-order transition as 
in Fig. El We note that in this case the holon density should be the same as the doublon 
density, because the holon and doublon densities must be the same to keep the density for 
the route of the bandwidth-control transition. 

To describe the both types of the transitions, we take the natural order parameter as 
(, where ( = X for the filling-control transition and ( = D for the bandwidth-control 
transition. Namely, we take 

F = J dr[\a (T - T c )C(rf + ^C(r) 4 - ^((r)]. (5) 

Here, fi^ is — U for the bandwidth-control transition whereas is the chemical potential \x 
conjugate to the doping concentration X for the filling-control transition. The Ising univer- 
sality is resulted from this Ginzburg-Landau- Wilson functional. When the Ising universal- 
ity is correct, we obtain the critical susceptibility exponent defined by Xc = ^F/d^ 2 }^ 1 oc 
(T — T c )~ 7 with 7=7/4 for two-dimensional systems, d = 2 and 7 ~ 1.24 for d = 3. The 
order parameter exponent defined below T c as (Q oc |T — T c \^ obtained from dF/d( = at 
)Uf = satisfies (3 = 1/8 and (3 ~ 0.325 for d = 2 and 3, respectively. The exponent with 
varying fi^ defined as (£) oc ^ 5 at T = T c is given by 5 = 15 and 4.8 for d = 2 and 3, 
respectively. These exponents for d = 3 were indeed claimed to be observed at the critical 
point of the bandwidth-control Mott transition for V2O3 26 . 

In the Ising-transition picture^, the transition is characterized by these simple exponents 
with the hyperscaling assumption being satisfied below the upper critical dimension d u = 4. 
The system has a single length scale £ which diverges at T = T c . In the present context, £ 
expresses the density correlation length or doublon density correlation length. 
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III. QUANTUM MOTT CRITIC ALITY 



A. General remark on quantum effect 

A nontrivial question arises when the critical temperature of the Mott transition, T c can 
be lowered. With lowering of T c , how do the quantum effects emerge? If T c becomes zero, 
then one might naively expect that conventional quantum critical phenomena would appear. 
A naive expectation would be that the transition might be described by the Ising universality 
class in d + 1 dimensions, where the additional one dimension emerges from the dimension 
in the imaginary time in the path integral formalism. It turns out later that this is not the 
case. In any case, the criteria for the existence of the non-negligible quantum effect should 
be determined from the existence of the Fermi degeneracy of the electrons. When the Fermi 
degeneracy temperature becomes comparable or higher than T c , the Mott transition has to 
be treated fully quantum mechanically. 

One may argue that even when the bare Fermi temperature is high, the effective Fermi 
temperature would be suppressed near the Mott critical point, because, at the continuous 
transition point to the insulator, the Fermi degeneracy temperature should be zero. However, 
the Fermi degeneracy may still coexist with the critical fluctuation in the metallic side near 
T c . In fact, even when T c is zero, the Fermi degeneracy always appears at temperatures 
sufficiently close to zero, if the parameter infinitesimally deviates from the critical point. 
Then the quantum effect should become relevant when one approaches the transition point 
by keeping the temperature sufficiently low. 

In the quantum region, we have to consider quantum dynamics. This can be done by 
considering the path-integral formalism, where the imaginary time direction must be addi- 
tionally considered in addition to the real spatial dimension. The time scale u)~ l diverges as 
uj^ 1 oc £ 2 in addition to the divergence of the spatial correlation length £. The quantum dy- 
namics is characterized by the dynamical exponent z. We note here that the Mott transition 
can be characterized by two different dynamical exponents in principle. This is because at 
the Mott transition, the single-particle spectra given by the single-particle Green's function 
and the two-particle correlations represented by the density (charge) correlation functions 
both have singular behaviors with diverging time scale, while these two may in principle 
follow different scalings. Therefore we can define two dynamical exponents z and z t for 
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single-particle and two-particle spectra, respectively. We will show below that these two 
coincide each other. 

Here, we discuss how the quantum effect alters the transition by assuming the region 
where T c is low or even zero. When T c becomes low, the Ginzburg-Landau expansion tells 
that the charge fluctuation becomes diverging accompanied by the quantum degeneracy, 
which becomes beyond the scope of the form (jSJ). One might expect that the quantum 
region could be described by the Ising universality in d + z dimensions. We will show that 
this does not apply and the quantum effect is more profound. 

B. What is different from the conventional quantum critical phenomena? 

In contrast to phase transitions with simple spontaneous symmetry breakings, the tran- 
sitions from metals to the band insulators and the Mott insulators have no spontaneous 
symmetry breaking by themselves. Therefore as in the gas-liquid transition, the metals and 
insulators have no clear distinction at nonzero temperatures, if the first-order transition 
would be absent. However, at zero temperature, insulators are always clearly distinguished 
from metals by the vanishing conductivity. Among various types of insulators, band insula- 
tors and Mott insulators both have clear distinction from metals by vanishing Drude weight 
and vanishing charge susceptibility (compressibility) at zero temperatur e 1 ' 34 . The Drude 
weight is the stiffness to the twist of the phase of the wavefunction in the spatial direction, 
while the charge susceptibility is the stiffness to the twist in the temporal direction in the 
path-integral formalis m 1 ' 32 ' 33 . Both of the two quantities have nonzero values only in metals. 
Therefore, metals under the perfectly periodic potential of ions are regarded as a state where 
symmetry of the phase of the spatially-extended electron wavefunction is broken. Then, at 
zero temperature, insulators cannot be adiabatically continued to a metal. Two phases have 
to be clearly separated by a phase boundary. 

From these facts, we notice that the Mott transition at zero temperature may have a 
quite different universality class. In fact, if one could lower T c in Eq. (J3J) by controlling some 
microscopic parameter, one may also expect that it could pass through zero and even to a 
negative temperature. This implies that the transition would become quantum critical and 
then the Mott transition would disappear as in the conventional scenario of the emergence 
of quantum critical phenomena as schematically illustrated in Fig. Ufa). However, we have 
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seen above that this cannot happen because of the clear distinction between metals and 
insulators at T = 0. This by itself indicates that the Ginzburg-Landau- Wilson scheme has 
to break down when T c becomes zero. When T c is lowered to zero, and if one tries to drive 
the control parameter further to lower T c , it in reality keeps T c at zero, namely the quantum 
transition at T = continues as we see in Fig. Efb). We call this continuation line, the 
T c = boundary. When T c becomes just zero from nonzero values, we call this point the 
marginal quantum critical point, which is indicated by the solid circle in Fig. Efb). 

We have discussed already that the first-order transitions with the critical end point are 
indeed found in experiments. We have also discussed and will discuss later that the metal- 
insulator transition through the T c = boundary appears to exist. Then the marginal 
quantum critical region is the crossing point of these two regions. We will clarify that 
many strongly correlated systems including the high-T c cuprates may be located in this 
marginal quantum critical region. Numerical results of the filling-control transition in the 
two-dimensional Hubbard model indeed suggest a continuous transition at zero temperature 
with the diverging charge susceptibility (density susceptibility) 21 i 22 i 2 ? , which is consistent 
with what is expected in the marginal quantum critical region. Therefore, this quantum 
criticality is certainly a realistic possibility. From experimental point of view, the quantum 
parameter g may be controlled by the lattice structure, particularly, by the geometrical 
frustration as already discussed in §1. The control of the frustration parameter was actually 
achieved by the choice of anions in K-type ET compounds while uniaxial pressure may also 
be used to control the frustration effects in general. 

C. Single-particle quantum dynamics 

The metallic phase except for one-dimensional systems has the adiabatic continuity with 
the Fermi liquid. Therefore, low-energy part of the free energy can be described by the 
fermionic operators of renormalized single particles, where the higher-order terms are renor- 
malized to the single-particle coefficient. In the insulating side as well, single-particle Green's 
function G describes the charge dynamics and may be given from a quasiparticle description 
with a gap A c (g) as 




(6) 
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FIG. 4: Schematic phase diagram of conventional quantum critical phenomena (a) and the Mott 
transition (b) in the parameter space of the control parameter conjugate to the order parameter, 
the parameter g to control the quantum fluctuation and temperature T. In (b), the T c = boundary 
continues as a line beyond the marginal quantum critical point. The shaded areas are the first-order 
phase boundaries. 

where the Matsubara frequency is u n , and E(q,k) = ±a/ A c (g) 2 + e(q, k) 2 with the bare 
dispersion e. In this expression, k is the momentum coordinate perpendicular to the locus 
of e(q, k) — and q denotes that parallel to e(q, k) = 0. Here E(q, k) is assumed to 
satisfy E(q, k) > 0(< 0) in the electron-doped (hole-doped) region (, namely for the pole 
of the upper (lower) Hubbard branch ). The imaginary part of the self-energy ImE and 
the renormalization factor Z is not considered here, because the singularities of the Mott 
transition are our main interest in this article while we assume that the singularities are not 
altered by ImE and Z^. 

Now we take the hole picture for the hole doping side so that E always takes E > both 
in electron and hole doped regions. Then, aside from the rigorous validity of the details 
of the above form for E, around the Mott gap edge, one can assume that the dispersion is 
expanded in terms of k as 

E(q,k) = a(q)k 2 + b(q)k 4 + (7) 

We have shifted the chemical potential to cancel the gap A c . Here, the gap edge is not 
necessarily isolated points in the momentum space, but may be a line or a surface, which 
may evolve to the Fermi surface in the Fermi liquid. The fc-linear term does not exist because 
we have assumed that k = is the gap edge: The fc-linear term violates the requirement 

14 



E > for small negative k. The cubic term should also be vanishing in the region of our 
interest because it becomes relevant only when a becomes sufficiently small, while then the 
cubic term also violates the initial assumption of E(q, k) > for negative k. The analyticity 
of the dispersion at small k is our assumption. This is actually plausible when the transfer 
energy for the distant pair of Wannier orbitals well converges to zero with increasing distance 
in the Hubbard-type models. 

The coefficients a and b are obtained as renormalized values after eliminating the higher 
order terms of the quasiparticle operator. In the metallic side, the rigid band picture is not 
justified. However, it is still legitimate to consider the quasiparticle dispersion around the 
Fermi level and the coefficients a and b as effective quasiparticle coefficients obtained in the 
evolution process of the metallic phase. This means that a and b may depend on the distance 
from the Mott transition point, whereas they still behave continuously. It should be noted 
that the variations of a and b can again be renormalized, which leads to ( independent a 
and b near the transition point. 

If the gap edge is given by isolated points and the coefficients a(q) at these gap edges 
approach nonzero positive constants on the verge of the transition, the dynamical exponent 
z characterized by the single-particle dispersion is given by z = 2 as in the generic transition 
to the band insulator— with a finite effective mass of quasiparticles. We see below that this 
does not hold any more at the marginal quantum critical point. 

From this quasiparticle description, the transition between metals and Mott insulators 
are described by the change in the quasiparticle dispersions. One of our central statements 
in this paper is that the criticality of the transition has one to one correspondence with the 
singularity of the quasiparticle dispersions. When the coefficient a stays positive through the 
transition at finite number of isolated points of gap edge, the free energy has a similar sin- 
gularity with the transition between the band insulator and metals, which is the continuous 
transition at zero temperature along the T = boundary in Fig^b). The first-order tran- 
sition may evolve only when a becomes zero. In the following, we clarify how the character 
of the transition and the single-particle dispersion are related each other. 
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D. Relation to Free Energy Form 



Here we relate the quasiparticle dynamics and the singularity of the free energy at the 
Mott transition in the quantum region. When we take the path integral formalism with the 
imaginary time r, the singular part of the free energy density is formally written by using 
the quasiparticle dispersion E as 

F = -(T/N)\nZ (8) 
Z = f n^(r)P0*(r)e- s /^, (9) 

rh/T rh/T 

S= / drV^^ 8 + / dr#(0*,0), (10) 
Jo i Jo 

where the effective hamiltonian H generating the single particle excitation E in Eq. ((7j) is 
rewritten by using the Grassmann variables <pi and </>* at the site i. 

The above quasiparticle form leads to the effective action at the chemical potential fi as 

S = ^2<j)*(q, k,u n )(-iu n - fji + E(q, k))(p(q, k,u n ) (11) 

i,n 

with the Matsubara frequency u n . 

Let us first study the filling-control transition. Although the expansion (jSJ) does not hold 
in the quantum region, the singular part of the free energy at zero temperature still has an 
expansion with respect to the doping concentration X. From the Matsubara- frequency and 
wavenumber dependent path integral form of the quantum dynamics leading from Eqs.(jHJ) 
to (jllj) . we obtain for the singular part of the free energy density at the transition, 

POO 

F = -Xfx-T dED(E)log(l + e- {E -^ /T ) (12) 
Jo 

with D(E) being the singular part of the density of states of quasiparticles. The Boltzmann 
constant is taken to be unity for our temperature scale. The particle density measured from 
the insulating phase is given by 



DC 



X= / dEf(E)D(E) (13) 
Jo 

with the Fermi distribution function f(E) = l/(e^ £;_M - ) / T + 1). 

We first consider zero temperature and the case where the dispersion E has minima at 
finite number of isolated points in the momentum space with the dispersion given by Eq.(|7|). 
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Then the particle density is given by 



X = A d k d F , (14) 

where A d is a dimensionality-dependent constant and kp is the Fermi wavenumber measured 
from the gap edge at the dispersion minima. Here we ignored the possible anisotropy of kp 
because it does not alter the essential part of the results below for the scaling properties. 
From the above relations, at T = 0, we have 

r-, , r 2a .X.2,1 46 X . 4 , i , 

= ~ + ^feW + ttJtJ 1 + • • •'• (15) 

We note that we are considering only the singular part of the free energy at the transition. 

When the dispersion proportional to a around these points qi in Eq.Q is present, the 
quasiparticle picture predicts that the total free energy measured from the insulator has the 
lowest order term at T = as 

F + Xfi oc aX^ d+2 ^ d . (16) 

From this free energy form, the charge susceptibility shows the scaling Xc = (d 2 F/dX 2 )~ 1 oc 
X 1_2 / d , which is the same as the transition to the band insulator. Then the first-order 
transition does not take place, because the charge susceptibility does not diverge for d > 2. 
The metal insulator transition is well defined only at zero temperature. This means that 
the transition occurs across the T c = boundary illustrated in Fig. 0] 

By starting from this continuous transition at zero temperature, the first-order transition 
can evolve in two fashions. One possibility is the case where a large Fermi surface satisfying 
the Luttinger theorem appears immediately upon doping. In this case, one has to take that 
Fermi surface as the locus E(q, k) — [A — and one gets 



by eliminating kp from 



and 



F + X/iocX 3 (17) 



"kp 

X oc I dk (18) 
o 



pk F 

F + X/ioc / k 2 dk, (19) 
Jo 

where we perform the integrations in the region around the locus E—fi = with the assump- 
tion a > everywhere. In this case, one gets the charge susceptibility Xc = {d 2 F / dX 2 )' 1 oc 
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X -1 , which means that the transition occurs at the marginal quantum critical point, where 
T c is still zero. This is because \ c diverges only at X = at zero temperature. To real- 
ize the phase separation (namely, the first-order transition), we need a further additional 
degeneracy of the dispersion at the gap edge. 

The second possibility is the marginal quantum critical point emerging with the vanishing 
a term at some isolated points go of the gap edge. In this case, we have the lowest order 
term 

F + Xfi oc bX {d+z)/d , (20) 

with z = 4, which yields 

Xc oc X 1 -^. (21) 

The exponent z = 4 appears because we are left with the quartic term proportional to b 
when the quadratic term proportional to a vanishes. 

E. Electron differentiation 

Now it turns out that the two possible ways of realizing the first-order transitions eventu- 
ally become merged to a unified picture, because even when we have the locus of E — \x = 
with Eq. (jI7J) being satisfied, a further flattening of the dispersion at the gap edge with van- 
ishing a is required to realize the first-order transition. It is unlikely that such a flattening 
emerges uniformly on the locus E = 0. Instead, it generically occurs from particular points 
of the locus because of the initial anisotropy of the Fermi surface in the band structure and 
the anisotropic correlation effects as well. Namely, the point with the smallest amplitude of 
a becomes zero first as a special point of the E — \i = surface. Then a quartic dispersion 
appears at this special point of the E — \i — surface when the system becomes marginal, 
which results in z = 4. 

Therefore, the Mott criticality of the marginally quantum critical point is also character- 
ized by an inevitable evolution of the electron differentiation, if the large Fermi surface is 
involved in the metallic side. The singular differentiation generates a quartic dispersion at 
particular points of the Fermi surface coexisting with dispersive generic part. 
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F. Comparisons with numerical and experimental results 



This large dynamical exponent z = 4 was suggested in several independent numerical 
calculations for the filling-control transition of the Hubbard and t-J models in two di- 
mensions at T = ^21. These are the exponent estimated from the compressibility in 
the form (|2T]) 21 i 22 i 23 i 37 . the Drude weight D estimated from the form D oc x 1+ ^ z ~ 2 ^ d ^^, 
single-particle dispersion 40 and the localization length £ oc (// — /ij -1 / 2 in the insulator side 
estimated from Green's function Eq.fJBJ), where £ is obtained from the Fourier-transformed 
spatially-dependent Green's function G(r,uj = 0) ~ exp(— r/^) 41 -. These imply that the 
Hubbard and the t-J models are located close to this marginal quantum critical point. 

Although the quasiparticle picture does not hold, it has numerically been shown that the 
dynamical exponent indeed becomes z = 4 at the marginal quantum critical point in the 
one-dimensional Hubbard model with next-nearest-neighbor transfers 4 ^ This suggests that 
the form (J2U|) is universally valid irrespective of the applicability of the quasiparticle picture. 

It is insightful to compare experimental results obtained for the high-T c cuprates with 
the present picture of the electron differentiation. In the high-T c cuprates, flat dispersions 
are universally observed near (tt, 0) and (0, tt) points in the angle-resolved photoemission 
experiments 43,44 . This flatness is beyond the conventional expectation obtained from the 
van Hove singularity, while it is a natural consequence, if the Mott critical point at T = T c 
is located at low temperatures. At and below T c , vanishing quadratic dispersion, given by 
a = should emerge in a region of the expected Fermi surface. 

The distance from the marginal quantum critical point in the phase diagram in Fig.0]may 
depend on details of materials and models. In addition, the actual go positions responsible 
for the marginal quantum Mott criticality may also depend on materials and models. For 
example, the change in relative amplitude of t' to t in the Hubbard model may change 
the location of go- This may even change the nature of the transition from the route across 
the T c = boundary to the route through the first-order transition illustrated in Fig. 0] This 
change may indeed occur from one high-T c compounds to another 24,45 . Actually the singular 
points go may deviate from (tt, 0) and (0, tt) points for larger t' . Width of the critical region 
may be influenced by the amplitude of t' as well. 

For the moment, definite assignments are not possible, but we infer two possibilities 
for the cuprates. One is that the system is indeed close to the marginal quantum critical 
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point and the Mott criticality is controlled by the flat dispersion. Even when the Mott 
criticality is controlled by these flat points, the experimental observation of the flat part of 
the dispersion at the Fermi level may have some difficulty because the strong damping is 
inevitably accompanied. The arc structure is observed in the angle-resolved photoemission 
experiments for the underdoped cuprates, which literally means that the Fermi surface 
around the flat-dispersion region is missing. This implies that the experimental resolution 
might not allow the detection of the Fermi surface around the flat part, because of the strong 
damping, while this flat part may govern the criticality. 

The other possibility is that the part of the flat dispersion is slightly away from the 
Mott gap edge. This is inferred from the fact that the flat-dispersion level in the hole- 
doped cuprates is slightly lower than the arc part around (tt/2, 7r/2), which is the dispersion 
minimum (in the hole picture) at the real Mott gap edge indicated by experiments^^ and by 
model calculations as well^. In this case, with the lowering doping concentration, the system 
first shows the marginal quantum critical behavior reflecting the flat dispersion. However, 
with further approaching the real critical point, it crossovers from the marginal quantum 
criticality characterized by z = 4 to the ordinary class z = 2, which eventually flows to the 
criticality for the T c = boundary. The real high-T c cuprates appears to have a variety 
between these two possibilities depending on the compounds. We will discuss consequences 
of the latter case further in Sec. IIIIKI 

Aside from these details and uncertainties, the overall structure of very different evo- 
lutions of the Fermi surface depending on the momenta is completely consistent with the 
picture that the transition metal compounds including the cuprates show electron differ- 
entiation arising from the proximity of the marginal quantum Mott criticality. Electron 
differentiation should become prominent when the system becomes closer to the marginal 
quantum critical point, while the electrons are more or less uniform along the T c = bound- 
ary. 

G. Breakdown of Ginzburg-Landau- Wilson scheme 

From the results in Sees. IIII CI and IIII D| the free energy near the marginal quantum 
critical point is generally expressed as 

F = -fiX + aX^l d + bX( d+ ^ d , (22) 
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where a and b have absorbed numerical constants in Eq.(|15|) as well as effects of the renor- 
malization factor Z being less than unity, and the constants have properly been rescaled. 

Even in the case of the bandwidth-control transition, when one can regard the closing 
of the gap by hole doping around a point g h an d simultaneous particle doping around qo p 
with the constraint of keeping the electron density n = 1, the above relation Eq. (J22j) may 
be replaced with D as 



because the "doublon" and "holon" concentrations are nothing but the above self-doping 
concentration of particles and holes. 

Now instead of Eq.(J5|). the free energy at zero temperature is expanded by ( and obtained 
after rescaling of the parameters as 



with the constraint ( > 0. It should be noted that this expansion of the free energy in 
terms of ( is obtained from the path integral form with the spatial as well as imaginary time 
dependence explicitly taken into account. This form of the free energy clearly violates the 
Ginzburg-Landau- Wilson scheme because the form of the free energy itself has d dependent 
nonanalytic expansion. 

The metal-insulator transition across the T c = boundary with a > is driven by the 
/i£ term in Eq. ()24j) . When fi^ is negative, the free-energy minimum exists at ( = 0, which 
corresponds to an insulator. Whereas the metallic phase is represented by the minimum at 
a nonzero positive (, which is realized by a positive The T c = boundary is determined 
from /i^ = 0. The criticality of this T c = boundary is given for the order parameter as 



In the "mean-field approximation", (3 is obtained from the spatially uniform derivative 
dF/dC = as 



F = UD + aD^' d + bD^/ d , 



(23) 



F = -/i C C + a( id+2)/d + bC^' d + cC {d ^ ),d ■ ■ ■ 



(24) 



(0 oc |/^| 



(25) 



(3 = d/2. 



(26) 



The susceptibility is given by 



XC = [^] 



-i 



2(d + 2) 
I 2 




(27) 
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leading to the susceptibility exponent denned by 

X< « |/i C |- 7 (28) 

in the metallic side /i( < as 

7 = 1- d/2. (29) 
At jj, — 0, we obtain another exponent from 

xc ~ or 1 ?- 5 (30) 

with 

5 = 2/d (31) 

in the metallic side. These are all "mean-field exponents" , although we will show these expo- 
nents are indeed correct. In the Ginzburg-Landau mean-field theory, the critical exponents 
do not depend on the dimensionality. However, this <i-dependent form of the free energy 
leads to the d dependent exponents even in the mean-field treatment. 

When a becomes zero, the critical point = becomes marginal and the first-order 
transition evolves if a becomes negative. This point with fj,^ = a = at T = is nothing but 
the marginal quantum critical point. Namely, the marginal quantum critical point at T = 
may be reached at a control parameter g = g c , for a = a^{g — g c (T)), > and b > 0. 

One might argue that the marginal quantum critical point looks similar to the conven- 
tional tricritical point^. because the continuous transition converts to the first-order transi- 
tion at T = at this point. However, it is qualitatively different because the Mott transition 
contains only the insulator and metal phases and no additional competitions a priori exist. 

Similarly to the transition across the T c = boundary, the critical exponents of the 
marginal quantum critical point have d dependent forms even in the mean-field treatment. 
The exponent f3 defined by the order parameter at g < g c and ^ = as 



is given by 



(C)<x\g-9c\ (32) 



f3 = d/2. (33) 



Near g = g c , \c, is expressed as 

X( = [g]-. „ [Hlf^fe _ 9c)< */- + ^pl b( ^ + HOt+V^-i. (34) 
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Then 

Xc^b-fcl -7 (35) 
holds in the metallic side g < g c , yielding the "mean-field" exponent 

7 = 2- d/2. (36) 

At g = g c , 

X( ~ b' 1 ?- 6 (37) 

with 

5 = A/d (38) 

is obtained. We note that the divergence of the susceptibility is stronger at the marginal 
quantum critical point than along the T c = boundary. It is also stronger for lower spatial 
dimensions. For example, 7 is 1 and 1/2 in two and three dimensions, respectively, for the 
marginal quantum criticality. On the contrary,7 is not positive in two and three dimensions 
along the T c = boundary. 

Even though the Ginzburg-Landau- Wilson scheme breaks down, the scaling relations 
(35 = 7 + (3 and a + 2(3 + 7 = 2 are satisfied both for the T c = boundary and the 
marginal quantum critical point. Along the T c = boundary we obtain a = 1 — d/2 while 
a = —d/2 for the marginal quantum criticality. They are consistent with the hyperscaling 
law 2 — a = (d + z t )v with v — 1/2 and zt = 4 for the marginal quantum criticality and 
z t — 2 along the T c = boundary. Here, z t should be the dynamical exponent of this Mott 
transition, namely the dynamical exponent for the density or doublon density correlations, 
which turns out to coincide with the dynamical exponent for the single particle excitations, 
z determined from the quasiparticle dispersion. 

We can also confirm that the hyperscaling relation is satisfied in the following way: The 
scaling relation and the exponents are derived from the scaling form of the free energy, 

F(a, M( )=r^/(ae 9 ,^) (39) 
with a scaling function / and the correlation length 

^(g-g c )- 1/2 ^C 1/d , (40) 

which implies y g = 2. Here the crossover exponent = 4 (y^ = 2) is derived from the 
dynamical exponent of the density fluctuations given by z t — 4 (z t = 2) for the marginal 
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quantum criticality (T c = boundary), respectively. The hyperscaling relation holds because 
this scaling form is satisfied. In fact this scaling form is derived from the single length 
scale £ which diverges at the transition point. This correlation length is indeed proportional 
to the mean carrier distance oc X~ l l d in the filling-control transition. This is obviously 
the single length scale which diverges at the transition. Satisfaction of the hyperscaling 
form (j23J) also clearly shows that this criticality is the consequence of the spatio-temporal 
quantum dynamics of two-particle excitations in the path integral form, where d and Zt 
represent spatial and imaginary-time fluctuations, respectively. 

H. Two-dimensional case 

In one and two dimensions, the powers of expansions in Eq. (j24j) stay at integers (for 
example, in ID, (d + 2)/d reduces to 3 and in 2D, (d + 2)/d reduces to 2). In two dimensions, 
the free energy is reduced to 



This again does not belong to the conventional scheme of the Ginzburg-Landau- Wilson for- 
malism, because the odd order term (the cubic term here) is not allowed in the conventional 
Landau expansion from the constraint of the symmetry around the critical point. This has 
some resemblance to the breakdown of the Ginzburg-Landau- Wilson scheme at the Lifshitz 
point of the structural transition although the physics contained here is quite different. Here 
the asymmetry is allowed because the part of negative ( does not exist. This can be easily 
understood in the analogy to the trivial metal to the band-insulator transition, where the 
carrier density cannot be negative either. Even in the transition between metals and band 
insulators in the noninteracting systems, the free energy has a similar form to Eq. (|24j) and 
the Ginzburg-Landau- Wilson scheme does not hold. 

In the Mott transition in two dimensions, the susceptibility is given by 



Then, 7 = 1,(3 = 1 and 5 = 2 hol d 24 ' 2 ^ . Remarkably, this agrees with recent experimental 
results on a k-ET compound, k-(ET) 2 Cu[N(CN) 2 ]C1 by Kagawa, Miyagawa and Kanod a 18l2 ° . 
The exponents of the finite-temperature critical point estimated by the conductance are 



F = -yUcC + ao(g - 9c)C 2 + K 3 + < 4 - 



(41) 




(42) 
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indeed consistent with these values for (3, 7 and 5 within the experimental accuracy. We 
note that this compound has a structure of highly two-dimensional anisotropy. For the 
moment, it is not well clarified how the crossover to the three dimensionality arising from 
weak three-dimensional coupling should appear experimentally. We will discuss below that 
these unusual exponents are also obtained practically even at finite temperatures, which is 
relevant in the realistic experimental condition of the finite-temperature Mott transition. 

I. Validity of the mean-field theory 

Very close to the marginal quantum critical point, the present mean- field exponent only 
marginally breaks down, because the Ginzburg criterion^ d+ zt > {2(3 + r y)/v with v = 1/2 
being the correlation length exponent indicates that the system is always at the upper critical 
dimension irrespective of the dimensionality, because the equality d + zt = (2[3+j)/v always 
holds. Here, v = 1/2 is a direct consequence of y g = 2. The fluctuation beyond the mean- 
field theory becomes irrelevant above the upper critical dimension. Although logarithmic 
corrections may exist, the mean- field description is thus basically correct at any dimension in 
this case. This explains why the "mean-field exponents" are observed in the k-ET compound. 
Below we restrict our analysis to the mean-field study because the primary exponents are 
correct. Although the mean- field exponents are correct, the hyperscaling also holds at any 
dimension. This peculiar compatibility is explained by the fact that the system is always at 
the upper critical dimension at any dimension. Detailed analysis of corrections based on the 
renormalization group study will be reported elsewhere. The same argument for the validity 
of the mean-field exponents is applied for the T c = boundary. 

Although the diverging density fluctuation is an inevitable consequence of the Mott crit- 
ical point for the marginal quantum criticality and the classical Ising criticality, it is highly 
nontrivial effect from the viewpoint of the weak coupling picture. In fact, naive perturba- 
tion expansions result in suppressions of the density fluctuations when the Mott transition 
is approached and the available perturbative treatment fails in reproducing this criticality. 
The one-loop calculation does not account for the Mott criticality, which in principle has to 
be explained in the real part of the self-energy to be calculated from the higher order loops, 
while it is so far an open issue to be derived in the future. In this sense, even the mean-field 
theory is not a straightforward framework in contrast with most of the mean-field theories 
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as those in magnetic transitions. 



J. Finite temperature effect 

The parameter g also has temperature dependence, because the quasiparticle dispersion 
in general has a temperature dependence. As we have clarified above, g is determined 
from the quasiparticle dispersions. Therefore g — g c can be replaced by T — T c , where the 
critical temperature becomes nonzero. This is the dominant finite-temperature effect at low 
temperatures. However, we have quite independent origin of the finite-temperature effects 
originating from the entropy term T( In £ in addition to the above temperature dependence in 
the free energy F. This generates an essentially singular contribution for /i^ < 0, because the 
free energy has the extremum at (o oc exp[/^/T]. This is an exponentially small contribution 
at low temperatures for ^ < 0. Therefore, it does not contribute to the present scaling 
behavior near T = in the power of T. At high temperatures, however, without the Fermi 
degeneracy, the expansion around the extremum reproduces the regular Ginzburg-Landau 
form 

A)(T-T C )C 2 , BC t 
Ft = - + — - McC, (43) 

when we redefine as ( — Co — ► C 

When T is nonzero, strictly speaking, the Ising universality class may appear in the 
region extremely close to the critical point even at low temperatures. However, this real 
critical region becomes exponentially narrow with decreasing temperatures and the quantum 
criticality governs out of this region. This crossover of the criticality is due to the essentially 
singular contribution of the entropy term. 

The results obtained for V2O3 26 indicates that T c is above 400 K and is high enough 
so that the Ising classical universality is well observed. On the other hand, T c for k- 
ET2Cu[N(CN)2]C1 is below 40K and low enough^, so that the quantum region dominates 
in the experimental results. All of these are consistent with what were observed. More 
quantitative analyses are left for future studies. 
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K. Effects of long-range part of Coulomb interaction for filling-control transitions 



The role of the long-range part of the Coulomb interaction between electrons is some- 
times controversial. In this subsection, we clarify how the long-range interaction modifies 
the conclusions obtained for the models only with the short-range interaction. Because the 
phase separation occurs for models with only the short-range force as the first-order transi- 
tion in the classical region, we restrict ourselves to the case of the Ising classical universality 
in this subsection, Sec JIIIKl Because of the long-range part and the resultant electrostatic 
condition, the real phase separation into two different electron densities is not allowed in the 
filling-control Mott transition. Reflecting the electrostatic condition, the diverging charge 
susceptibility at strict zero wavenumber is eventually suppressed for the filling control tran- 
sitions. However, it still causes critical fluctuations at nonzero and small wavenumbers. 

The Poisson equations for the external test charge p ext and the induced charge p in d are 
given by 

7 2 i Pext 



J ext j 

e 



V 2 md = -?f, (44) 

where <p ext is the electrostatic potential generated by the external test charge, while the 
induced charge generates additional potential <pi n d- The dielectric constant is e. When the 
external test charge p ex t is at the origin with the unit charge, in metals, the Thomas-Fermi 
screening occurs as p inc i = —a<p to t for weakly /c-dependent term with the definition of the 
total electrostatic potential (f) t ot — 'Pext + 'Pind- in fact, by the Thomas-Fermi screening, a is 
estimated as a = 2meqF/(irh 2 ), where m, e and qF are the electron effective mass, charge and 
the Fermi wavenumber, respectively. Here we have assumed isotropic sphere of the Fermi 
surface for simplicity. Since the external test charge at the origin leads to a wavenumber 
independent form p e xt{q) = 1 in the momentum representation, we obtain 

1 

<Ptot = — 2 (45) 

a + eq 2 

and the resultant induced charge as 

—a 

Pind = — ■ o, (46) 

a + eq 2 

which indicates the Yukawa-type normal screening, pi n d{f) = aexp[— 2y/irear}/r in the real 
space. When one considers the susceptibility as the linear response induced by the imposed 
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external charge as the perturbation, the susceptibility x defined by /w = —xPext is now- 
given by 



X = -— r (47) 



a 

a + eq 2 

Whereas, if the Coulomb potential is regarded as the perturbation, the linear response 
defined by p ind = ~Xc<\>ext is given by 

X, = (48) 
a + eq 1 

which shows that the charge susceptibility is suppressed near q = because of the long-range 
Coulomb force. With this screening, the free energy in the classical region is modified to 

F(r) = J dr'^Xc 1 ^ ~ r')X(r') 2 + h)X(rf - ecpX, (49) 

where the Fourier transformed susceptibility is given by 

X:\q) = X cH 1 (?) + — 2 , (50) 
eq 2 

= - + c(a)q 2 . (51) 
a 

Here, XcH is the charge susceptibility for the Hubbard model without the long-range Coulomb 
part. The second term in the right hand side of Eq. ()5()j1 represents the long-range Coulomb 
interaction. The scaling of a near the critical point is determined from the effective mass 
m and is given as a oc £ 4-d and the function c(a) scales as c(a) oc £ 2 ~ d . The charge 
susceptibility is enhanced except for extremely small wavenumber. 

In the realistic condition, the minimum of Eq.(|49jl appears at q m in ~ l/(ec) 1//4 . In the 
cuprates, this is roughly estimated to be q m i n ~ 0.17T, since e is the order of 10-100, and 
c ~ 10 as we see below. Namely, the charge susceptibility becomes strongly enhanced in 
the region of small wavenumbers, which are in the order of one tenth of the Brillouin-zone 
size, 2tt. Even around the marginal quantum Mott criticality, this induces a "softening" 
of the charge response with poor screening in the nm length scale, which causes strong 
dynamical fluctuations of electron density at these small wavenumbers. This provides us 
with mechanisms of various unusual properties for metals near the Mott insulator. If F 
becomes negative at q m in, an instability toward the charge ordering occurs. 

Even when T is below T c in the Hubbard-type models, where only the short-range force 
is considered, the phase separation dynamics toward the q = mode in reality freezes at a 
stage of a finite q constrained from the long-range repulsion for the filling-control transition. 
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It should be noted that in the bandwidth-control transition, the divergence at q = 
indeed occurs. This is because the electrostatic condition is not violated for the diverging 
doublon susceptibility. 

We now realize that the larger dielectric constant induces the larger enhancement of 
the charge susceptibility as obtained by substituting q min into Eq. (|5L)p . Actually, if good 
metallic carriers independently exist in addition to the strongly correlated electrons which 
yields the doped Mott insulator, such good metallic carriers efficiently screen the long-range 
Coulomb interaction of the correlated electrons. This allows more possibility for the part of 
the correlated electrons to approach the marginal quantum critical point. In fact, when the 
good metallic carriers perfectly screen the motion of the correlated electrons, it corresponds 
to the limit e — ► oo in the above argument and the diverging density susceptibility for the 
correlated part of electrons can indeed occur at q — 0, because the density fluctuations 
are completely compensated by the good-metallic carriers except for the onsite interaction. 
This may occur in a two-band system, where one band has large bandwidth supplying good 
metallic carriers, and the electrons on the other band are strongly correlated near the Mott 
insulating phase. Although it is not a simple two-band system, the situation in the hole- 
doped cuprate superconductors is in a sense ideal from this viewpoint, because upon the 
carrier doping, "good metallic carrier" first appears around the (vr/2, 7r/2) region of the Fermi 
surface, while the carriers near (it, 0) region with strong correlation effects become doped in 
the presence of these itinerant carriers. Then the instability for the phase separation or the 
inhomogeneity actually occurs as the inhomogeneity of the u (tt, 0) carriers" compensated by 
the "(7r/2,7r/2) carriers". To realize this situation, it appears to be important to recognize 
that the cuprates are located close to the marginal quantum critical point, but strictly 
speaking located slightly in the side of the T c = boundary. 

Even when the charge ordering or the phase separation does not occur, the spatial inhomo- 
geneity is easily driven by impurity potential or lattice distortions, because of the underlying 
enhanced density susceptibility. This may be relevant as the mechanism of the structure 
observed by scanning tunnel microscopy (STM) in the cuprates and mang 
The spatial inhomogeneity in the long length scale (typically at 1 to 10 nm scale) observed 
experimentally cannot be explained by the naive Thomas- Fermi screening length Atf, since 
nominally Atf is the order of 1/qF ~ 0.3nm. The inhomogeneous structrure may also be 
enhanced by the experimental condition probed at the surface. The present results obtained 
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from the Mott criticality has a tight connection to the approach from dynamical stripe fluc- 
tuations^^, while the importance of the underlying Mott criticality has not been recognized 
in the literature. 

IV. NON-FERMI-LIQUID PROPERTIES 
A. Mode coupling theory 

Now we discuss the consequences of the enhanced density susceptibility near the marginal 
quantum critical point. Within the mean- field theory, we assume the dominant part of the 
susceptibility in the vicinity of a small and nonzero momentum Q and around the zero 
frequency as 

= -fa, + D,(IO + fa -Qf + ■■■) ■ (K) 
where D s is the diffusion constant of the density fluctuations, K is the distance from the 
marginal quantum critical point, and T is a constant. Near the critical point, they follow 
the critical scaling as oc K oc £ _1 and D s oc £~ 2 , which reproduces the above 
scaling Eq. (j35|) with Eq. (j40j) at uj = 0. This again satisfies the dynamical exponent z t = 4 
because u scales as D S K 2 oc £ -4 , while the scaling of K reproduces v = 1/2. Since it is at 
the upper critical dimension, this Ornstein-Zernike-type form is justified. We note that the 
enhancement is much stronger for d = 2 than d = 3. 

It should be noted that the characteristic energy scale of this charge fluctuation is much 
larger than that of spin and orbital fluctuations, because c(a)q 2 should typically have the 
energy scale of the Mott gap at the boundary of the Brillouin zone. This is clear because at 
the zone boundary, the density fluctuation requires generation of the spatially alternating 
doublons and holons. Then D s has the energy scale comparable to the Mott gap. For the 
filling control transition, it may have an energy scale as large as several eV as in the case of 
the cuprates (~ 2 eV), while for the bandwidth-control transition, this energy scale is not 
necessarily large. The large energy scale of these Mott fluctuations explains many puzzling 
properties of metals near the Mott insulator as we will discuss below. 

Since the dominant fluctuation occurs at small but nonzero wavenumbers, the conserva- 
tion law as in the density conservation at Q = does not exist and the dynamical exponent 
stays at Zt — 4 for the filling-control transition. For the bandwidth-control transition, the 
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FIG. 5: Temperature dependence of the inverse of the density susceptibility peak. A "Curie- 
Weiss" temperature dependence is observed in an extended region for several choices of doping 
X = 0.1,0.2 and 0.3. This linear temperature dependence causes non-Fermi-liquid properties. 

real divergence at q = again does not involve the conservation law because the doublon 
density does not commute with the Hamiltonian and is not conserved. Therefore, the dy- 
namical exponent again stays at four and the form Eq. (j52j) is justified. This is in contrast 
with the ferromagnetic fluctuation at Q = for the spin fluctuations, where iu in Eq. (|52|) 
is replaced with iu/'E with 5 oc q so that the dynamical exponent increases by oneS. 

Now we formulate a mode coupling scheme for the electron-density and doublon-density 
fluctuations originating from the Mott criticality around small Q. The inverse of the static 
susceptibility is obtained from Eq. (j4^|) as 

xdQ^ 1 = 2a (T - T c ) + 66(Cq) + 12c ^(C 9 2 ), (53) 

q 

where we introduced the average of the density fluctuation as ((g). This is given from the 
fluctuation dissipation theorem as 

1 f°° 

(Q = 7T / ducoth(u/ZT)lm X dq,u). (54) 

In the inverse of the bare susceptibility, the mass term is renormalized by the mode- 
coupling term proportional to Y^qiCq) m Eq. (j53|) . By renormalizing the zero temperature 
value, T c is renormalized to 

K = T C --J2(Q(T = 0), (55) 
a - 
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where 



(C g 2 )(T = 0) = - / dulmxdq^) 



1 f°° 



(56) 



is the zero-temperature fluctuation. Then the renormalized inverse susceptibility is given by 



We can solve Eqs. (|52|) . (J57)l and (J58)l selfconsistently, where we take into account the 
Gaussian fluctuations of the density through the mode coupling. In the mode-coupling 
scheme, the critical exponents stay at the mean-field values, which is justified from the 
above arguments in Sec. 1111 11 However, effects of the fluctuations are taken into account in 
a selfconsistent fashion. A similar phenomenological theory has been formulated by Moriya 
for spin fluctuations 55 . 

B. Perturbative self-energy 

In the following, we employ a perturbative treatment to understand the non-Fermi-liquid 
behavior as well as superconducting instability. The perturbative scheme has a limited 
applicability for these highly nonperturbative phenomena of the density fluctuations. How- 
ever, when the Mott criticality is properly taken into account through the phenomenological 
treatment of the density fluctuations, we expect that the essence can be captured even when 
we employ a perturbation theory for the other part of calculations. 

In the actual calculations for the realistic choice of parameters, a reasonable set of pa- 
rameters can be derived in the following: We choose the parameter values of a filling-control 
transition for a two-dimensional system appropriate for the copper oxide superconductors 
inferred from the frequency dependence of the optical conduct ivit y 56 > 57 , characteristic size of 
the observed inhomogeneit y 48 1 49 and the doping dependence of the density susceptibility in 
numerical^ and experimental results, which suggest a ~ 0, b ~ 0.7, c ~ 100, T^ 1 ~ 3X 
and D s ~ 30X by taking the energy unit t(~ 0.4eV) and the lattice constant as the 
length unit. The characteristic wavenumber and energy of the density fluctuations are 
roughly 7r/10 and 0.5-leV, respectively, which determines b and D s , respectively. In fact, 




(57) 



with the definition for the finite-temperature correction; 




(58) 
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the dielectric function e(q = 0,u) obtained from the optical conductivity provides us with 
Imx(g = 0,o') = Im[l/e(g = 0,oo)] — 1. The obtained results for Imx deduced from the 
experimental data for the optical conductivity have a prominent peak structure around 0.5- 
1 eV 57 , which indicates the characteristic energy scale of the density fluctuations. After 
considering the screening effect by the long-range Coulomb part, this suggests D s has the 
order of leV. The uncertainty of the parameters remains because of the lack of accurate 
experimental probes to estimate the frequency and wavenumber dependences of dynamical 
density fluctuations. Basically, all the results presented here do not depend on Q within the 
choice < Q < 0.2tt. 

As an example, we consider the carrier doping in the Hubbard model with the dispersion of 
the square lattice E(q) = — 2t(cos q x +cos q y ) with additional input of the density fluctuations 
given by Eq. (j52j) . Through the mode coupling, the solution of the selfconsistent equations 
. (|57 |) and (|58|l shows that the Curie- Weiss type behavior X( ~ (T + 0)" 1 holds in an 
extended temperature region with small Weiss temperature near the marginal quantum 
critical point. FigureElshows calculated results of such Curie- Weiss bahaviors in an extended 
temperature region. The criticality stays at the mean-field form, while the Curie- Weiss form 
is retained over very large temperature region with renormalized values of coefficients. Even 
when ao = is employed, we obtain the linear temperature dependence of namely the 
Curie- Weiss form in a wide temperature region because of the linear temperature dependence 

of HQ- 

Now we calculate the electron self-energy. The electron self-energy in the perturbation 
expansion up to the second order of the interaction for the filling control is given as 

TU 2 

E(q, u n ) = ^2 G ( k > iuJ n)xx(q - k, i(u n - u rn )). (59) 

k,n 

Here the imaginary part of the self-energy ImE is governed by the Curie- Weiss behavior of 
Xx through Eq.(|o^jl. Because of the linear temperature dependence in x^ 1 in an extended 
temperature region, we obtain the linear temperature dependence also for ImE in a wide 
temperature region. In the marginal-critical region with xx oc X~ x , the standard Fermi- 
liquid behavior ImE oc T 2 is replaced with the non-Fermi-liquid form ImE oc T, which 
may cause various unusual properties. The resistivity in two-dimensional systems becomes 
nearly proportional to T as p oc ImE oc T in contrast to the standard Fermi-liquid scaling 
p oc ImE oc T 2 . 
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A long-standing puzzle in the doped Mott insulators is widely observed long-tail structures 
in the optical conductivity extending up to the order of 1 eV in various transition metal 
oxides and organic conductors ^ 56 ' 57 ! 58 . The tail has a structure of power law decay in the 
optical conductivity as <j(uj) ~ oj~ p with p ranging between 0.3 and 1. Origin of such long 
tail structure has to be attributed to fluctuations in the energy scale of 1 eV and cannot be 
accounted for by the spin and orbital fluctuations, since they have much lower energy scale 
typically less than 0.1 eV. We note that the density fluctuation mechanism examined in this 
article naturally accounts for such fluctuations at large energy scale. 

V. SUPERCONDUCTIVITY EMERGING FROM MOTT CRITIC ALITY 

A. Pairing originating from the marginal quantum Mott criticality 

It is widely recognized that the origin of the high-temperature superconductivity in the 
copper oxides® has to be explained by considering the strong electron correlation effects, 
although the mechanism is still puzzling and not definitely figured out. After the discovery, 
various aspects of magnetic mechanisms were extensively examined. From the weak coupling 
pictur o 59 ! 60 , the spin fluctuation theories were considered, where strong antiferromagnetic 
fluctuations were assumed to mediate the Cooper pairing. Then the origin of the high- 
T c superconductivity was assumed to arise from the criticality of the antiferromagnetic 
quantum critical point. In the strong-coupling expansion represented by the t-J model 6 ^, 
it was claimed that the pairing is basically through the singlet formation stabilized by the 
superexchange term proportional to J in the t-J model. In both of the approaches, the 
mechanism of the Cooper pairing is more or less the same and they are categorized as the 
magnetic mechanism. 

On the other hand, from the initial stage of the studies on the cuprate superconductors, it 
has been well recognized that the superconductivity occurs in the region of the doped Mott 
insulator near the Mott transition 62 . However, since the criticality of the Mott transition 
was not well identified until recently, the role of the Mott criticality for the mechanism 
of the superconductivity was not well appreciated. In fact, as we already clarified in the 
previous sections, the Mott transition itself is a transition driven by the order parameter of 
the electron (or doublon-holon) density, and has nothing to do with the magnetic degrees 
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of freedom by itself. Although the antiferromagnetic fluctuations occur at (tt, it) in the 
magnetic Brillouin zone for the square lattice, the Mott criticality occurs independently 
through quite a different fluctuation, namely through the singular density fluctuations at 
small wavenumber. If the antiferromagnetic order exists at low temperature of the Mott 
insulating phase, the antiferromagnetic fluctuations around (it, tt) coexist with the density 
fluctuations around small wavenumber arising from the Mott criticality. It should also be 
noted that this fluctuation has a completely different origin from the ordinary charge-order 
fluctuations at a commensurate wave vector, although the present instability may also trigger 
the charge ordering. 

It is naturally expected that the density fluctuations at small wavenumber inherent to 
the Mott transition may play a novel role in stabilizing the superconducting phase. In fact, 
the density fluctuations may be the origin of instabilities to various symmetry breakings 
including not only superconductivity but also charge and magnetic orderings, since the 
diverging density fluctuations are directly connected with the flattening of the quasiparticle 
dispersions at the Fermi level leading to the enhanced density of states. The diverging 
density of states widely enhances the instability for various orders. 

The pairing mechanism arising from this density fluctuation was recently studied by 
a perturbative scheme of the mode-coupling theory^. We here discuss superconductivity 
assuming the proximity to the Mott quantum critical point in detail. In the present scope, 
the density fluctuations have the largest energy scale over spin and orbital fluctuations and 
are the primary origin of the unusual properties. This proximity of the Mott transition 
indeed mediates the Cooper pairing through the enhanced density fluctuations. When we 
follow the perturbative scheme, the effective interaction between two electrons is obtained 
from the density fluctuations as 



up to the second order in U with xx obtained from ([52)1 . Of course other fluctuations as 
spin fluctuations also affect the effective interaction. However, to extract the role of density 
fluctuations clearly, we ignore the contribution from spin fluctuations in this article as the 
first step. In fact, the density fluctuation plays the dominant role because of its large energy 
scale. When spin fluctuations are also considered, we expect that it reinforces the pairing, 
since both of the fluctuations enhance the same type of pairing symmetry while they do not 




(60) 
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FIG. 6: Eigenvalues A for the linearized Eliashberg equation (|61|) are plotted as functions of 
temperature. Since the eigenvalue for the d x 2_ y 2 symmetry is the largest, it is plotted for X = 0.1. 
The transition temperature T sc is obtained from the temperature where the eigenvalue crosses 
A = 1 (horizontal line). We have used the density susceptibility in Fig. [5J We plot for several 
choices of the Hubbard U interactions with the energy unit of the nearest-neighbor transfer t on 
the square lattice Hubbard model. 

interfere each other because of their fluctuations at very different wave numbers as we see 
below. We obtain the linearized Eliashberg equation for the superconducting gap A as 



where iV is the number of sites. This Eliashbrg equation is solved selfconsistently. Con- 
sidering the level of perturbative approximations here, the first nontrivial way to solve this 
problem is to take the bare Green's function for G, and we ignore the normal self-energy 
corrections to G. In the calculation of Green's function in Eq. (j61|) . the standard Hubbard 
model on the square lattice is employed as an example. However, the nonperturbative effect 
is taken into account through Eq. (j52j) for xx with the parameter values introduced above. 
Then Eq. fl6()j) is inserted to Eq. (|61j) and the eigenvalue A is calculated for the right-hand 
side of Eq. (j61|) . Namely, the linearized Eliashberg equation is solved selfconsistently for the 
relevant parameter values for the cuprate superconductors as cited above. 




(61) 
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B. Unconventional pairing 

The solution of the Eliashberg equation Eq.(|61|) for the parameter values above shows that 
the right hand side of the linearized Eliashberg equation has the largest eigenvalue for the 
d x 2_ y 2 pairing symmetry, which leads to the highest superconducting transition temperature 
T sc for this pairing symmetry as we show in Fig. |H1 Figure [7| shows that the d x 2_ y 2 pairing 
symmetry indeed wins over the eigenvalues for the other symmetries including the symmetry 
of the extended s-wave-symmetry pairing. It is remarkable that even though the fluctuations 
are at a small wavenumber Q, it generates the anisotropic pairing. This is because the 
effective interaction is repulsive in the most part of the Brillouin zone because of the first term 
in Eq. (|6U|) . while it becomes attractive only in the small wavenumber region. This forces the 
pairing to have an anisotropy with nodes. The d x 2_ y 2 symmetry and its node position are 
understood because the largest gap grows in the (ir, 0) and (0, 7r) regions, which is stabilized 
by the flat dispersion. Then the only possibility is to make nodes in the diagonal direction 
in the Brillouin zone. We note that the gap amplitude may be substantially underestimated 
because we have underestimated the flatness of the dispersion in Green's function by taking 
the bare Green's function instead of the correct one. Figure |B] shows how the eigenvalue 
grows with lowering temperatures. The superconducting transition temperature within this 
approximation is estimated from the temperature where the eigenvalue exceeds unity in 
Fig. |H1 The transition temperature has the order of O.Olt to 0.05t as we see in Fig. which 
corresponds to the order of 100K for the copper oxides when we take t ~ 0.4eV. It should 
be noted that the large energy scale (namely the Mott gap scale ~ 2 eV) of the density 
fluctuation represented here by the parameter value D s = 30X is crucial for achieving such 
a high transition temperature. In other words, the high-energy excitations substantially 
contribute to enhancing the transition temperature. It is remarkable that within this simple 
approximation, the quantum Mott criticality has a dramatic effect on the superconductivity, 
which is comparable or even larger than that by the magnetic mechanism in the same level 
of approximations. 

If we properly consider the self-energy effects in Eq. (|61|) . we expect that the supercon- 
ducting transition temperature T sc becomes vanishing at the Mott transition point although 
the pairing interaction is most enhanced at the Mott transition point because of the critical 
enhancement of the density fluctuations. These two should cause the separation of T sc and 
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FIG. 7: Temperature dependence of the eigenvalue for the linearized Eliashberg equation for the 
d x 2_ y 2 (squares) and extended s-wave (circles) symmetries. 
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FIG. 8: U dependence of the superconducting transition temperature for the d x 2_ y 2 symmetry in 
the energy scale of t. The transfer is assumed to have t ~ 0.4 eV and for other parameter values 
see the text. 

the gap amplitude. It leads to the pseudogap behavior in the underdoped region. This is left 
for future studies. We note that one has to be careful in taking account of the self-energy 
effects because a part of it appears through the density fluctuation itself, which is already 
taken into account here but is beyond the presently available perturbative treatment in the 
literature. 
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The present analyses based on the perturbative treatment and mean-field-type Eliashberg 
equation do not take into account fluctuation effects particularly in two-dimensional systems. 
In fact, in purely two-dimensional systems, we expect Berezinskii-Kosterlitz-Thouless (BKT) 
type transition for the gauge symmetry breaking of the superconductivity and this aspect 
is not considered here. Nonetheless, the present analyses have significance in the following 
points: First, even when the BKT transition is expected in pure two-dimensional system, 
its transition temperature has a comparable value to the mean-field results as known in the 
analyses of the XY model universality. Therefore, the present mechanism of the supercon- 
ducting transition may also work for the BKT transition at the similar temperature scale, 
which can be inferred from the present simple approximations. Second, the significance 
of the present approach is that the superconducting mechanism arises from a completely 
new origin of the proximity from the quantum Mott criticality with the enhanced density 
fluctuations at small wavenumber. As a first step, clarification of possible relevance in the 
experimental situation is desired even at the mean-field level. Third, it is useful to com- 
pare consequences of the present mechanism with the conventional ones including magnetic 
fluctuation mechanisms at the same level of approximations. For the mean-field analyses, 
it is possible because we have many available results for the conventional mechanisms in 
the literature. Fourth, in the cuprates, very weak interlayer coupling, which still does not 
destroy the dominance of two-dimensional Mott criticality in a wide region, may sensitively 
induce real superconducting transition as the three-dimensional one. This circumstance may 
show that the mean-field treatment for the superconducting transition offers a qualitatively 
correct way of understanding if the pairing mechanism is correctly picked up. More detailed 
analyses with consideration of the inherent two-dimensional fluctuations and the resultant 
BKT transition is left for future studies. 

The present pairing mechanism may also work near the bandwidth-control transition 
point. In fact the basic mechanism can be straightforwardly applied to the region near the 
critical end point of the bandwidth-control transition. This may explain the superconducting 
phase observed near the Mott transition point of k-ET compounds family^. 

The superconductivity near the valence instability point was studied theoretically as a 
model for CeCu 2 Si 2 , CeCu 2 Ge 2 , and other heavy fermion compounds^^. Since we expect 
a similar novel criticality to the present Mott criticality, it would be an intriguing issue to 
pursue the mechanism of the superconductivity along the same line. 
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Since electron density fluctuations of course strongly couple with phonon modes, the 
phonons also supplement this density fluctuation mechanism. In fact, the strong coupling to 
small wavenumber phonon necessarily occurs, if the enhanced electron density fluctuations 
already exist in the Mott critical region. This may reinforce the strong coupling to phonons 
even when the conventional electron-phonon coupling A is rather weak. In fact, the B lg 
out-of-plane phonon mode with a small momentum transfer— may have a relevance and 
presumable resultant kink structure in the angle-resolved photoemission spectra should be 
considered under this circumstance. 

VI. SUMMARY AND DISCUSSION 

In this paper, it has been shown that the Mott transition is successfully described by a 
new framework for quantum phase transitions. The natural order parameter for the Mott 
transition is the electron doping concentration for the filling-control transition. For the 
bandwidth-control transition, the natural order parameter is the doublon (or holon) density. 
At zero temperature, the phase boundaries of metals and Mott insulators always exist as 
the Mott transition, which occurs either as the first-order or continuous transitions. 

If the Mott transition occurs as the continuous transition at zero temperature and other 
spontaneous symmetry breakings are not involved, the metals and insulators are adiabati- 
cally connected at finite temperatures. We call this regime the quantum regime surrounding 
the T c = boundary. On the other hand, if the transition occurs as the first-order transition 
at T = and terminates as the critical point above the Fermi degeneracy temperature, we 
call it the classical regime. Sandwiched by these two regimes, the marginal quantum critical 
region appears, where the first-order boundary and the continuous T c = boundary meet 
at T = 0. 

Although the density fluctuation is completely suppressed in the Mott insulator because 
of the Mott gap, the criticality of the continuous Mott transition in the metallic side can be 
described by the critical enhancement of the density fluctuation at small wavenumbers in con- 
trast to the naive expectation. This critical enhancement indeed occurs at finite-temperature 
Mott critical line as well as at the marginal quantum critical point. The marginal quantum 
criticality shows nontrivial and novel features. 

In the classical regime around the high critical temperature, the Mott criticality is de- 
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scribed by the Ising universality class. However, in the quantum regime, the clarified Mott 
criticality indicates the breakdown of the Ginzburg-Landau- Wilson (GLW) scheme. The 
quantum dynamics is derived from the path-integral formalism with the one-to-one corre- 
spondence between single-particle dynamics and two-particle correlations. The divergence 
of time scale at the transition is described by the dynamical exponent z = 2 for the quan- 
tum region and z = 4 for the marginally quantum region. Remarkably, it is shown that 
the free energy has nonanalytic expansion with respect to the order parameter, with power 
depending on spatial dimensionality d for the quantum as well as the marginally quantum 
criticality, in marked contrast with the GLW expansion, which is demensionality indepen- 
dent. This unusual expansion results in dimensionality-dependent critical exponents, which 
also indicates that the Mott transition occurs always at the upper critical dimension at any 
d. Then, at any d, the scaling relations and the hyperscaling are still satisfied while the 
mean-field description is basically justified for the critical exponents except for the possible 
logarithmic corrections. They are totally described by a new universality class. Particularly 
when the Mott critical temperature becomes lowered just to zero temperature, the marginal 
quantum critical point appears and the critical exponents are given by 7 = 2 — d/2, (5 = d/2 
and S = A/d. 

The present theoretical framework for the quantum Mott criticality has clarified many 
aspects which are consistent with the experimental results near the Mott insulator. These 
are summarized in the following: 

First, the universality class of the quantum Mott transition well explains the otherwise 
puzzling critical exponents recently discovered in k-(ET) 2 Cu[N(CN) 2 ]C1; 7 = = 1 and 
5 = 2. This exponents are identified as those at the marginal quantum critical point in 
two-dimensional systems. The scaling description is completely consistent with both the 
classical Ising- type transition observed in V 2 3 and the quantum transition observed in 
k-(ET) 2 Cu[N(CN) 2 ]C1 as well as in /t-(ET) 2 Cu 2 (CN) 3 . 

Second, such unusual exponents and the new universality class at the marginal quan- 
tum critical point inevitably cause the differentiation of electrons in the momentum space 
even when the large Fermi surface with the contained Luttinger volume is expected. This 
draws a concrete picture how the Fermi liquid breaks down to the Mott insulator. The 
differentiation along the Fermi surface is the driving mechanism of the emergence of the flat 
dispersion, and the arc structure observed by the angle resolved photoemission experiments 
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in the cuprate superconductors. The differentiation generates particular points on the Fermi 
surface responsible for the criticality, which is the reason why the hyperscaling relation is 
satisfied in the present theory. 

Third, approaching the marginal quantum critical point, the system becomes more and 
more sensitive to the external perturbations and is easily driven to the inhomogeneous state, 
which has been suggested in various types of surface probes. We have estimated the typical 
length scale of the inhomogeneity determined from the balance of the Mott criticality and 
the long-range Coulomb interaction in the filling-control transition. 

At the same time, the universality class of the Mott criticality at finite temperatures 
in the classical region is protected from the randomness because the Ising universality at 
finite temperatures is not influenced by the small randomness. The first-order transition 
is only driven by the mechanism of the Mott transition equivalent to the Ising class. The 
universality class of the critical point at the termination point of the first-order transition 
may receive fewer effects of Anderson localization, because the Anderson transition does not 
drive the first-order transition at all. The insensitivity to the randomness is particularly true 
in three dimensions. In purely two dimensional systems, however, the situation is nontrivial 
because of the sensitive random field effects on the Ising transition, which will be discussed 
in a separate paper. This is also the origin of the tendency for the spatial inhomogeneity. At 
the marginal quantum critical point, we also expect that the universality is protected against 
randomness at least in three dimensions because of diverging density of states at the gap edge 
and efficient screening through divergent density susceptibility. Near the marginal quantum 
critical point, the enhanced density of states implies that the Thomas- Fermi screening length 
becomes short, which leads to a more efficient screening of the random potential by the fewer 
carriers. Then contrary to the naive expectation, the density fluctuation at a long length 
scale appears. Along the T c = boundary, the metallic phase may be under a severe effects 
of randomness and the continuous metal-insulator transition is eventually triggered by the 
Anderson localization. 

The filling-control transition requires a special care because of the long-range Coulomb 
interaction. For charged electrons, the real divergence of the density fluctuation at strictly 
q = is suppressed. However, the density fluctuations are still strongly enhanced at a small 
nonzero wavenumber, which may cause unusual properties. In addition, if dispersive light 
carriers coexist with the carriers near the Mott insulator, the screening and compensation 
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by such dispersive and good metallic carriers allows the critical divergence of the density 
fluctuation for the correlated carriers. This appears to be realized indeed in the cuprate 
superconductors. 

The fourth point for the experimental relevance of the marginal quantum criticality is the 
non-Fermi-liquid properties. The non-Fermi-liquid properties observed in the doped Mott 
insulators with two-dimensional anisotropy are accounted for by the density fluctuations 
arising from the quantum Mott criticality. This in fact explains the T-linear resistivity in 
two dimensions, when we employ the mode coupling theory. In particular, the fluctuations 
at the energy scale as large as 1 eV observed as the long tail in the optical conductivity 
universal in transition metal compounds and organic compounds are explained not by spin 
or orbital fluctuations but by this density fluctuation of the quantum Mott criticality. 

We have also shown that the mechanism of the high-T c superconductivity can be ascribed 
to the density fluctuations originating from the marginal quantum Mott criticality. Although 
the density fluctuation occurs at small wave number, it causes the unconventional pairing 
with a nodal structure. For the realistic parameter values of the cuprate superconductors, 
the solution of the linearized Eliashberg equation has the highest transition temperature 
of the order of 100K for the d x i_ y i wave symmetry . It is remarkable that, within the 
present level of approximation, the density fluctuations inherent near the quantum Mott 
criticality overlooked in the literature cause comparable or even larger effects than the spin 
fluctuation mechanisms extensively studied for the cuprate superconductors. The large 
energy scale of the density (charge) fluctuations may help even stronger instability towards 
the superconductivity, if this instability could be more carefully tuned by the design of 
material parameters. This is a challenging future task. 

From the experimental point of view, it is highly desired to develop a good experimental 
probe for studying the dynamical and short-range density (charge) correlations. In contrast 
to the magnetic correlations well studied by neutron scattering and NMR, experimental 
probes for the wavenumber and frequency dependent charge correlations of electrons are 
poor. We can study the optical conductivity and dielectric functions only at zero wavenum- 
ber while STM and several microscopes can detect only the static structure on the surface. 
Raman scattering is not powerful enough so far for the study of the systematic wavenumber 
dependence. Electron energy loss spectra and inelastic X-ray scattering in principle probe 
the dynamical density correlations while the present energy resolution is rather poor. To 
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uncover the whole feature of physics near the Mott transition, it would be highly desired 
to develop experimental probes for the frequency and momentum dependence of electron 
density correlations and it will make a breakthrough in this field. Our prediction from the 
present work is that the extended charge (density) fluctuations at small wavenumbers may 
be observed in the frequency and wavenumber dependent spectra in the critical region of a 
certain class of the Mott transitions in transition metal compounds and organic conductors. 
The density fluctuation may also occur with a compensation between two different types of 
carriers, good metallic carriers and the strongly correlated carriers near the Mott insulator. 
This compensated density fluctuation will require a more refined probe to be detected. This 
will reveal the missing ring of various puzzling properties including the mechanism of the 
high-T c cuprate superconductors and the criticality of the Mott transition in the organic 
conductors. 
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